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The color glass condensate approach describes successfully heavy ion collisions at RHIC. We investigate Iron-air 
collisions within this approach and compare results to event generators commonly used in air shower simulations. 
We estimate uncertainties in the extrapolation to GZK energies and discuss implications for air shower simulations. 



1. Introduction 

Higher twist corrections to hadronic interac- 
tions become increasingly important at high en- 
ergies. An effective way to resum all these con- 
tributions is the color glass condensate picture. 
We review the Kharzeev-Levin-Nardi (KLN)[T] 
approach to heavy ion collisions, which has suc- 
cessfully been applied to RHIC physics, and ex- 
trapolate this model to Iron-air collisions at ener- 
gies up to the GZK (w 10 19 7 eV) cutoff. So far, 
we only calculate multiplicities. In the near fu- 
ture, we hope to present a full model which also 
treats the forward scattering in detail, which is 
important for air showers. 




2. Review of BBL 1.0 

In Ref. [2] we introduced the black body limit 
(BBL) model for hadron-nucleus reactions. Va- 
lence quarks scatter coherently off the gluon field 
in the target. The transverse momentum ac- 
quired in this reaction is of the order of the satura- 
tion momentum, defined by the density of gluons 
in the target. This leads to independent frag- 
mentation of the leading quarks at high energies, 
when the saturation momentum is high. The 
leading baryon effect known from lower energies 
is suppressed. 

Gluon production in this approach is realized 
within the KLN model: a simple ansatz for the 
unintegrated gluon distribution function (uGDF) 
is applied to the fc t -factorization formula. In this 
paper, we apply this model to Iron-air collisions. 



Figure 1. The construction of the uGDF for low 
densities. The full line shows the uGDF of a single 
nucleon. The dashed line is the average uGDF for 
Pa ~ 0.5. The dotted line shows the uGDF for 
the averaged low density. 



3. KLN approach to nucleus-nucleus colli- 
sions 

In the fc^-factorization approach [3J, the distri- 
bution of produced gluons is given by 

dN g _ AN C <Pp± f p± , 2h 

<Pr ± dy ~ N*-lJ pi J dk±as 

x cj) A (x 1 ,k 2 ± )<j) B (x2,{P±- k ±) 2 ) (!) 
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Figure 2. The multiplicity of charged particles 
as a function of the centrality. The grey dotted 
line shows the result when taking Q 2 S to be pro- 
portional to the average density Ta- The data is 
from PHOBOS (5). 



with N c = 3 the number of colors. Here, p± 
and y denote the transverse momentum and the 
rapidity of the produced gluons, respectively. 
The light-cone momentum fractions of the col- 
liding gluon ladders are then given by £1,2 = 
p± exp(±y)/y / s, where y/s denotes the center of 
mass energy and y is the rapidity of the produced 
gluon. We set p™ ax such that the minimal satu- 
ration scale QJ lln (xi i 2) in the above integration is 
K Q cd = 0.2 GeV. 

The KLN approach [T] employs the following 
uGDF: 

1 Q 2 

(f>(x,k x ;r±) j-^- s 2 , (2) 

a s {Qi) max(Qj,k ± ) 

where Q s denotes the saturation momentum at 
the given momentum fraction x and transverse 
position rj_. The overall normalization is deter- 
mined by the multiplicity at mid-rapidity for the 
most central collisions. The saturation scale for 
nucleus A is taken to be proportional to the den- 
sity of participants, 7ip art (rjJ. This is not a uni- 
versal quantity which depends only on the proper- 
ties of a single nucleus, in other words, the uGDF 
is not fully factorizable. 



Figure 3. The saturation scale as a function of the 
rapidity y — log(l/a;). The vertical lines denote 
the forward and mid-rapidities for the different 
energy regions. 



A possible solution seems to be to define the 
saturation momentum squared to be proportional 
to Ta{ti_), which is dependent on one nucleus 
only and therefore respects factorization. How- 
ever, this ansatz cannot describe the data on the 
multiplicity as a function of centrality for Au-Au 
collisions at RHIC, as shown in Fig. [2] Why this 
does not work can best be seen with the following 
example. We consider a peripheral Au-Au colli- 
sion, and want to construct the uGDF at the edge 
of one of the nuclei. Here, Ta is very small, which 
means that only in some nucleus configurations, 
we actually find a nucleon at this position. Let 
us denote the probability to find (at least) one 
nucleon with pa- The uGDF at this position is 
then pa times the uGDF of a single nucleon, as 
sketched in Fig. [1] (dashed line). Taking the av- 
erage density Ta would lead to the dotted line, 
and gives the wrong uGDF, since averaging has 
to be done after constructing the wave function 
and not before. The density under the condition 
to find at least one nucleon is Ta/pa, see [4] for 
details. The saturation scale is therefore: 
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Figure 4. The rapidity dependence of the multi- 
plicity for central (b — 2.4 fm) Pb-Pb collisions 
at the LHC. 



and the ansatz for the uGDF is 
c(> A =p A <p(^) . (4) 

In Ref. [4] we explain that the multiplicity is a 
homogeneous function of order one in the density 
of both nuclei, and the n par t description is a good 
approximation of the factorized KLN approach. 

4. Results 

Before showing some results and extrapolating 
to high energies we want do discuss the evolu- 
tion of the saturation scale as a function of the 
energy. Fig.[3]shows Q s as a function of the rapid- 
ity for fixed coupling and running coupling evolu- 
tion. See Ref. [2] for details on these two types. 
The two evolution types show large differences 
at forward rapidity for LHC and GZK energies. 
Differences at mid-rapidity for RHIC and LHC 
are however less significant. Therefore, fixed and 
running coupling predict very similar results for 
RHIC and LHC overall multiplicities. 

4.1. Accelerator Results 

Very good agreement with the data from the 
PHOBOS collaboration [5] is shown for the multi- 
plicity at mid-rapidity as a function of centrality, 



Figure 5. The multiplicity of charged particles as 
a function of the centrality for a 5500 GeV Pb- 
Pb collision at the LHC (b=2.4 fm). The upper 
curve for each A shows fixed coupling, the lower 
curve shows running coupling evolution. 



as shown in Fig. [2j As already stated, the differ- 
ences between fixed and running coupling evolu- 
tion are small. 

In Fig. [4] we see predictions for central 
(b=2.4 fm) Pb-Pb collisions at the LHC (^/i = 
5500 GeV) . The centrality dependence of the mul- 
tiplicity at mid-rapidity is shown in Fig. [5] 

4.2. Iron Air 

Fig.[6]shows the rapidity dependence of charged 
particles for central Iron-Nitrogen collisions at the 
three reference energies for RHC,LHC and GZK 
(200 GeV, 5500 GeV and 400 TeV, respectively). 
Up to LHC energies, differences at mid-rapidity 
are smaller than 15%. Only at GZK energies, we 
see a qualitative difference. 

Fig. [7] shows the total multiplicity of charged 
particles as a function of the lab-energy per nu- 
cleoli. Since in cosmic ray physics the energy 
of a nucleus is usually the total energy, plotting 
Eiab/A is a useful quantity when comparing to 
proton-air collisions. We also compare to the 
standard hadronic interactions models Sibyll 2.1 
[3] and QGSJET-IIc [7]. First, we observe that 
nucleus air collisions in Sibyll are obtained by su- 




Figure 6. The rapidity dependence of charged 
particles for central Fe-N collisions at the three 
reference energies. The results for QGSJET-II 
and KLN (fixed and running coupling evolution) 
are shown. 



perposition of hadron air collisions. The multi- 
plicity scales therefore with the number of partic- 
ipants in the projectile and ranges from a factor 
30 to 36 compared to proton air collisions. In 
the QGSJET-II model and the KLN approach, 
screening in the initial state reduces the multi- 
plicity. The predicted results of these two ap- 
proaches are quite similar up to LHC energies, 
whereas above this energy, QGSJET-II predicts 
higher multiplicity than the color glass approach 
(running coupling evolution). 

4.3. Conclusions 

We compared the multiplicities of Iron-air col- 
lisions in the color glass condensate approach 
(KLN model) with hadronic interaction models 
used in air shower simulations. The results of 
QGSJET-II and KLN are quite similar up to LHC 
energies. Above this energy, QGSJET-II predicts 
higher multiplicities. 
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